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A  modeling  study  of  benzene  and  phenyl  radical  formation  was  performed  for  three  low-pressure  pre¬ 
mixed  laminar  flat  flames  having  an  unsaturated  C2  or  C3  hydrocarbon  fuel  (acetylene,  ethylene,  and 
propene).  Predictions  using  three  published  detailed  elementary-step  chemical  kinetics  mechanisms  were 
tested  against  molecular  beam  mass  spectrometry  (MBMS)  species  profile  data  for  all  three  flames.  The 
differences  between  the  predictive  capabilities  of  the  three  mechanisms  were  explored,  with  an  emphasis 
on  benzene  formation  pathways.  A  new  chemical  kinetics  mechanism  was  created  combining  features  of 
all  three  published  mechanisms.  Included  in  the  mechanism  were  several  novel  benzene  formation  reac¬ 
tions  involving  combinations  of  radicals  such  as  C2 1 1  +  C4H5,  C2H3  +  C4H3,  and  C5H3  +  CH3.  Reactions 
forming  fulvene  (a  benzene  isomer)  were  included,  such  as  C3H3(propargyl)  +  C3H5(allyl),  as  well  as 
fulvene  to  benzene  reactions. 

Predictions  using  the  new  mechanism  showed  virtually  all  of  the  benzene  and  phenyl  radical  to  be 
formed  by  reactions  of  either  C3H3  +  C3H3  or  C3H3  +  C3H5,  with  the  relative  importance  being  strongly 
dependent  on  the  fuel.  C5H3  +  CH3  played  a  minor  role  in  fulvene  formation  in  the  acetylene  flame.  The 
(  ’2 1 1 ,  +  C4Hx  reactions  did  not  contribute  noticeably  to  benzene  or  phenyl  radical  formation  in  these 
flames,  sometimes  being  a  major  decomposition  channel  for  either  fulvene  or  phenyl  radical.  The  formation 
pathways  for  C3H3  and  C3H5  were  delineated  for  the  three  flames;  although  the  key  reactions  differed 
from  flame  to  flame,  1CH2  +  C2H2  <=>  C3H3  +  H  was  important  for  all  three  flames. 


Introduction 

A  molecular-level  understanding  of  combustion 
chemistry  is  essential  for  controlling  the  products  of 
molecular-weight  growth,  which  can  lead  to  the  for¬ 
mation  of  polycyclic  aromatic  hydrocarbons  (PAHs), 
many  of  which  are  known  mutagens  and/or  carcin¬ 
ogens.  The  PAHs  themselves  are  likely  precursors  of 
soot  particles  in  flames.  For  flames  of  aliphatic  fuels, 
all  of  the  aromatic  rings,  including  the  first,  need  to 
be  formed  by  combination  of  smaller  hydrocarbon 
fragments,  thus  providing  motivation  for  under¬ 
standing  benzene  formation  in  combustion. 

Considering  the  level  of  interest  in  benzene  for¬ 
mation,  the  lack  of  consensus  on  what  the  major  re¬ 
action  pathways  are  might  be  surprising  to  the  casual 
observer.  Thorough  reviews  of  all  but  the  most  re¬ 
cent  work  have  been  provided  by  Westmoreland  et 
al.  [1]  and  Melius  et  al.  [2],  Two  broad  types  of  ring¬ 
forming  reactions  appear:  (1)  acetylene  (C2H2)  ad¬ 
dition  to  a  C4  vinyl  radical,  with  subsequent  cycli- 
zation,  that  is,  reactions  Rf  [3]  and  R2  [4]  in  Table 
1  and  (2)  combination  of  resonantly  stabilized  free 
radicals  (RSFRs),  with  subsequent  rearrangement 
and  ring  formation,  for  example,  reaction  R3a,  pro- 
pargyl  (H2CCCH)  radical  recombination,  2C3H3  <=> 
C6H5(phenyl)  +  II  [2,5,6]. 


In  the  acetylene-addition  reactions,  one  reactant 
(C2H2)  is  a  major  product  of  fuel-rich  combustion, 
present  in  high  concentrations,  and  the  adduct  does 
not  require  intramolecular  rearrangements  prior  to 
cyclization.  Plowever,  both  n-C4H3  (HCCHCCH) 
and  n-C4H5  (CH2CHCHCH)  have  more  stable 
isomers — i-C4H3  (H2CCCCH)  and  i-C4H5 
(CH2CHCCH2) — into  which  they  can  readily  be 
converted  [5],  greatly  diminishing  the  amount  of  the 
n-C4H3  and  n-C4H5  available.  While  RI  and  R2  are 
strongly  exothermic,  they  also  have  large  decreases 
in  entropy,  becoming  less  favored  at  higher  tem¬ 
peratures.  Furthermore,  R2  has  a  net  decrease  in 
number  of  moles,  causing  it  to  be  less  favored  at 
lower  pressures. 

Relatively  high  concentrations  of  C3H3  were  also 
seen  in  flames,  due  to  the  resonance-enhanced  sta¬ 
bility  of  propargyl.  In  low-pressure  acetylene  [1], 
1,3-butadiene  (C4H6)  [4],  and  benzene  [7]  flames, 
peak  mole  fractions  of  C3H3  near  10“ 3  were  ob¬ 
served.  Recombination  of  propargyl  radical  also  has 
no  barrier  to  forming  a  ring-containing  C6  species, 
whether  it  be  fulvene,  benzene,  or  phenyl  +  H 
[2,5,6];  R.3  can  therefore  become  a  direct  aromatic 
ring  formation  route,  especially  at  low  pressures. 
However,  reaction  R3a  has  a  crossover  temperature 
(above  which  AG  >  0)  of  1729  K,  and  is  much  less 
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ELEMENTARY  REACTION  KINETICS 


TABLE  1 

Rate  coefficients  for  reactions  which  are  discussed  in  the  text  and  in  subsequent  figures,  and  which  are  used  in  the 
expanded  mechanism  called  Mia  later  in  the  text 


Number 

Reaction 

A 

b 

E 

Reference 

Rl 

C2H2  +  n-C4H5  <=>  C6H6  +  H 

8.209  X  103 

0.801 

6348. 

[26],  pw 

R2 

U-C4H3  +  C2H2  <=>  C6H5 

1.667  X  1()10 

0.446 

7719.3 

[26],  pw 

R3a 

2  C3H3  <=>  C6H5  +  H 

5.00  X  1012 

0 

0 

pw 

R3b 

2  C3H3  <=>  fulvene 

4.50  X  1012 

0 

0 

pw 

R3c 

2  C3H3  <=>  n-C5H3  +  CH3 

5.00  X  1011 

0 

0 

pw 

R4 

C3H3  +  C3H5  <=>  fulvene  +  H  +  H 

5.562  X  102° 

-2.535 

1692. 

[9] 

R5a 

i-C5H3  +  CH3  <=>  fulvene 

1.40  X  1013 

0 

0 

[11],  pw 

R5b 

;-c5h3  +  ch3  <=>  c6h5  +  h 

6.00  X  1012 

0 

0 

[11],  pw 

R6a 

n-C5H3  +  CH3  <=>  fulvene 

1.00  X  1013 

0 

0 

[11],  pw 

R6b 

n- C5H3  +  CH3  <=>  C6H5  +  H 

1.00  X  1013 

0 

0 

[11],  pw 

R7 

+  C2H2  <=>  C3H3  +  H 

1.80  X  1014 

0 

0 

[5,21] 

R8a 

C2H3  +  CH3  <=>  C3H6 

3.688  X  1030 

-5.034 

7691.9 

[31] 

R8b 

C2H3  +  CH3  <=»  C3H5  +  H 

2.861  X  1021 

-2.118 

11841.8 

[31] 

R9a 

n- C4H3  +  C2H3  <=>  fulvene 

4.00  X  1012 

0 

0 

pw 

R9b 

n-C4H3  +  C2H3  <=>  C6H5  +  H 

1.60  X  1013 

0 

0 

pw 

RlOa 

n- C4H5  4-  C2H  <=>  fulvene 

4.00  X  1012 

0 

0 

pw 

RlOb 

n- C4H5  +  C2H  <=>  C6H5  +  H 

1.60  X  1013 

0 

0 

pw 

Rlla 

i-C4H3  +  C2H3  <=>  fulvene 

1.00  X  1013 

0 

0 

pw 

Rllb 

i-C4H3  +  C2H3  <=>  C6H5  +  H 

6.00  X  1012 

0 

0 

pw 

Rile 

i-C4H3  +  C2H3  <=>  2C3H3 

4.00  X  1012 

0 

0 

pw 

R12a 

i-C4H5  +  C2H  <=>  fulvene 

1.00  X  1013 

0 

0 

pw 

R12b 

i-C4H5  +  C2H  <=>  C6H5  +  H 

6.00  X  1012 

0 

0 

pw 

R12e 

i-C4H5  +  C2H  <=>  2C3H3 

4.00  X  1012 

0 

0 

pw 

Units  are  cm3,  mol,  s,  kcal.  pw  =  present  work. 


thermodynamically  favored  than  Rl,  with  a  cross¬ 
over  temperature  of  2218  K. 

Other  pathways  to  benzene  have  been  proposed. 
Lindstedt  and  Skevis  [8]  considered  C2H3  + 
C4H4(CH2CHCCH)  <=>  CgHg(benzene)  +  H  and 
C2H3  +  C4H6(CH2CHCHCH2)  <=>  C6H8(l,4-cy- 
clohexadiene)  +  H,  with  the  C6H8  decomposing  to 
benzene  and  H2.  Marinov  et  al.  [9]  proposed  a  com¬ 
bination  of  two  C3  RSFRs,  propargyl  and  allyl  (R4). 
Moskaleva  et  al.  [10]  considered  methyl  +  cyclo- 
pentadienyl  (also  a  RSFR),  forming  methylcyclopen- 
tadiene,  which  forms  fulvene,  then  benzene.  Mebel 
et  al.  [11]  calculated  the  potential  energy  surface  of 
C6H6(benzene)  <=>  C5H3  +  CH3  reactions.  The  re¬ 
verse  reactions  can  also  be  possible  benzene  for¬ 
mation  pathways  (R5,  R6)  from  either  of  the  two 
most  stable  forms  of  C5H3,  H2CCCCCH  (here 
called  i-CsH3)  and  HCCCHCCH  (n-C5H3).  Walch 
[12]  calculated  the  potential  energy  surface  of  both 
R2  and  the  reaction  of  i-C4H3  with  acetylene. 

There  are  not  infinitely  many  benzene  formation 
reactions.  Considering  only  bimolecular  reactions, 
the  two  reactants  should  collectively  contain  pre¬ 
cisely  six  carbon  atoms  and  roughly  six  hydrogen  at¬ 
oms  (for  the  above  examples,  five  to  nine  hydrogen 
atoms).  An  ultimate  goal  would  be  to  determine 


which  of  the  conceivable  benzene  formation  path¬ 
ways  are  important.  The  answer  would  be  complex, 
probably  depending  on  the  fuel  and  on  the  reactor 
conditions,  even  for  the  rather  restricted  case  of  low- 
pressure  premixed  flames  of  aliphatic  fuels. 


Approach 

To  assess  the  possible  importance  of  various  ben¬ 
zene  formation  pathways,  we  performed  a  modeling 
study  on  three  sets  of  flame  data.  The  flames  se¬ 
lected  were  all  fuel-rich  low-pressure  premixed  lam¬ 
inar  flat  flames  with  feeds  containing  an  unsaturated 
C2  or  C3  hydrocarbon  as  a  fuel,  oxygen,  and  argon; 
all  three  flames  had  species  concentration  profiles 
(including  benzene)  measured  by  molecular  beam 
mass  spectrometry  (MBMS),  as  well  as  experimental 
temperature  profiles.  The  three  flames  were  (FI)  a 
20  Torr  C2FI2  flame  with  45  mol  %  Ar,  cold  gas  ve¬ 
locity  of  97  cm/s,  and  an  equivalence  ratio  ( <f> )  of  2.5 
[13];  (F2)  a  20  Torr  C2H4  flame,  50  mol  %  Ar,  62.5 
cm/s,  <j>  =  1.9  [14];  (F3)  a  37.5  Torr  C3H6  (propene) 
flame,  25  mol  %  Ar,  48.2  cm/s,  <f)  =  2.32  [15],  Mod¬ 
eling  calculations  employed  the  CHEMKIN  pro- 
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grams  [16-20],  Full  multicomponent  diffusion,  in¬ 
cluding  thermal  diffusion,  was  calculated.  The 
catalytic  boundary  condition  for  H-atom  recombi¬ 
nation  at  the  burner  surface  was  also  included.  Ex¬ 
perimental  temperature  profiles  were  used  as  input, 
in  lieu  of  solving  the  energy  equation. 

Three  different  chemical  kinetics  mechanisms 
were  tested.  The  first  (Ml)  was  the  mechanism  of 
Pauwels  et  al.  [21],  derived  from  Miller  and  Melius 
[5],  plus  47  reactions  describing  the  chemistry  of 
propene  and  of  the  three  C3H5  isomers,  from  Mar- 
inov  et  al.  [22,23],  which  allowed  modeling  of  the 
propene  flame  (F3).  The  second  mechanism  (M2) 
was  from  Richter  et  al.  [24],  based  on  Zhang  and 
McKinnon  [25],  with  extensive  PAH  formation 
chemistry  at  low  pressure  (20  Torr).  To  afford  a  more 
complete  comparison,  the  reaction  1CH2  +  C2H2 
<=>  C3H3  +  H  (R7)  was  added  to  M2.  The  third 
mechanism  (M3)  was  from  Marinov  et  al.  [9,22,23]. 
M3  built  on  the  Miller  and  Melius  [5]  chemistry  by 
including  additional  reactions  of  RSFRs  and  PAH 
formation  reactions  at  atmospheric  pressure.  While 
many  other  chemical  kinetics  mechanisms  are  avail¬ 
able  in  the  literature,  the  three  chosen  here  provided 
a  reasonable  basis  for  comparison,  showing  areas  in 
which  Ml  can  be  improved. 

Modeling  Results 

All  three  mechanisms  did  well  at  predicting  the 
decay  of  fuel  and  02  in  the  three  flames.  Table  2  lists 
the  experimental  and  the  predicted  peak  mole  frac¬ 
tions  for  the  three  flames,  as  well  as  the  peak  loca¬ 
tions  (in  centimeter  height  above  burner).  The  pre¬ 
dictions  also  included  those  from  an  extensively 
modified  and  expanded  version  of  Ml,  called  Mia. 
The  key  changes  made  in  forming  Mia,  and  the 
modeling  results  using  Mia,  are  discussed  below. 

All  three  mechanisms  predicted  both  the  peak 
concentration  and  the  shape  of  the  C2H2  profiles 
quite  well  in  all  three  flames.  M3  overpredicted  peak 
C6H6  mole  fractions  for  all  three  flames.  Many  of 
the  pressure-dependent  reactions  in  M3  had  rate  co¬ 
efficients  that  are  fit  for  atmospheric  pressure,  so  M3 
overpredicted  molecular-weight  growth  in  such  re¬ 
actions  as  C2H3  +  CH3  <=>  C3H6  (R8a),  hence  the 
overprediction  of  C4H3  and  C4H5  in  all  three  flames. 
However,  C@H5  concentrations  were  underpre¬ 
dicted. 

Peak  CgHg  mole  fraction  predictions  from  M2 
were  within  25%  for  FI,  a  factor  of  10  low  for  F2, 
and  roughly  100  times  too  low  for  F3.  M2  did  not 
include  the  propargyl  recombination  reaction  (R3), 
which  most  likely  explains  the  CgHg  shortfall  and 
C3H3  overprediction  for  the  three  flames.  In  FI,  for¬ 
mation  of  C6H6  was  by  C3H3  +  C3H4(allene)  <=> 
CgHg  +  H,  yet  C3H4  was  also  overpredicted. 

The  peak  predicted  C6H6  mole  fractions  from  M2 
were  74%  of  the  peak  measured  value  for  FI,  11% 


for  F2,  and  0.89%  for  F3.  M2  does  not  include  the 
propargyl  recombination  reaction  (R3),  which  most 
likely  explains  the  C6H6  shortfall  and  C3H3  over-pre¬ 
diction  for  the  three  flames.  In  FI,  formation  of 
CgHg  was  by  C3H3  +  C3H4  (allene)  <=>  C6H6  +  H, 
yet  C3H4  was  also  overpredicted. 

The  predicted  peak  C0H6  mole  fractions  using  Ml 
ranges  from  30%  to  50%  of  the  data.  Peak  C3H3 
mole  fractions  were  overpredicted  for  all  three 
flames,  although  for  FI  and  F2  the  Ml  predictions 
came  closest  to  the  data.  Peak  C4H3  mole  fractions 
were  higher  than  the  data  for  all  three  flames.  Pre¬ 
dicted  C4H5  peaks  were  3.3%  of  the  measured  peak 
for  FI,  830%  for  F2,  and  84%  for  F3,  C6H5  predic¬ 
tions  were  12%  of  the  measured  peak  for  F2  and 
140%  for  F3,  following  the  same  trend  as  the  M2 
predictions. 


Expanded  Benzene  Formation  Chemistry 

Building  on  the  mechanism  Ml  [5,21],  a  signifi¬ 
cantly  expanded  mechanism  (Mia)  was  created,  con¬ 
taining  91  species  and  527  reactions.  (Note:  The 
mechanism  is  available  by  contacting  the  authors  at 
cjpope@sandia.gov  or  jamille@ca.sandia.gov.)  Be¬ 
yond  the  C3H5  and  C3H6  species  and  reactions,  two 
other  large  additions  were  made  using  reactions 
from  either  M2  [24]  or  M3  [22,23],  Enhanced  ben¬ 
zene  oxidation  chemistry  led  to  inclusion  of  five  C5 
species  (C5H5O,  C5H6,  linear  C5HS,  C5H40,  and 
C5H4OH).  Molecular- weight  growth  chemistry  up  to 
C10H8  (naphthalene)  was  added;  although  many  of 
these  reactions  were  contained  in  both  M2  and  M3, 
rates  from  M2  had  been  fitted  for  low  pressures. 

In  addition  to  the  sweeping  changes  mentioned 
above,  several  individual  reactions  were  added,  and 
several  key  rate  coefficients  were  modified.  Rate  co¬ 
efficients  for  reactions  R1  and  R2  were  from  Wang 
and  Frenklach  [26],  As  an  upper  bound  for  these 
rates,  the  20  Torr  fits  for  all  four  product  channels 
[26]  were  summed  and  used  as  the  rate  for  R2.  Sim¬ 
ilarly,  the  sum  of  the  four  product  channels  for  C2H2 
+  n  =  C4H5  for  0-90  Torr  [26]  was  used  for  the 
R1  rate  coefficient.  The  resulting  sums  were  fit  into 
the  three-parameter  Arrhenius  form.  The  barrier 
height  used  [26]  for  R2  was  3.1  kcal/mol;  the  barrier 
calculated  by  Walch  [12]  was  8.3  kcal/mol,  close  to 
the  AMI  barrier  height  [26],  Accordingly,  the  acti¬ 
vation  energy  for  R2  was  increased  by  5.2  kcal/mol, 
as  was  the  fitted  activation  energy  for  Rl.  The  rate 
for  C6H5  +  C2H2  <=>  C6H5C2H  +  H  was  taken  from 
Yu  et  al.  [27]  by  fitting  the  sum  of  the  rates  for  all 
product  channels. 

The  overall  rate  for  C2H3  +  C2H3  was  changed 
to  9.0  X  1013  cm3/mol  s  [28,29],  Consistent  with  the 
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TABLE  2 

Experimental  peak  mole  fractions  for  selected  species,  with  predicted  values  from  mechanisms.  Ml,  Mia,  M2,  M3: 

(a)  flame  FI,  (b)  flame  F2,  (c)  flame  F3 


(a)  C0H0 

flame — FI 

Species 

Ml 

Mia 

M2 

M3 

Data 

C2H2 

(0.55) 

(0.57) 

(0.60) 

(0.58) 

(0.56) 

c2h3 

4.78E-05(0.70) 

4.66E-05(0.80) 

1.74E-05(0.73) 

6.66E-05(0.75) 

3.77E-05(0.58) 

c3h3 

5.63E-04(0.72) 

5.28E-04(0.72) 

1.29E-02(0.82) 

1.89E-03(0.72) 

2.24E-04(0.80) 

c3h5 

1.27E-0.5(0.69) 

7.21E-06(0.59) 

1.70E-04(0.56) 

4.34E-05(0.71) 

xxx 

c4h3 

9.45E-05(0.80) 

9.20E-05(0.80) 

3.51E-04(0.88) 

1.02E-04(0.82) 

8.76E-06(0.59) 

C4H5 

1.04E-06(0.80) 

5. 77E -06(0.78) 

1.60E-06(0.69) 

6.94E-05(0.82) 

3.13E-0.5(0.79) 

c6hs 

2.24E-06(0.60) 

2.71E-06(0.62) 

5.23E-07(0.80) 

1.44E-05(0.69) 

xxx 

c6h6 

2.31E-0.5(0.55) 

4.99E-05(0.65) 

4.35E-05(0.59) 

7.27E-04(0.60) 

5.87E-0.5(0.60) 

(b)  CoH4 

flame — F2 

Species 

Ml 

Mia 

M2 

M3 

xData 

c2h2 

2.91E-02(0.95) 

3.03E-02(1.00) 

2.56E-02(0.75) 

3.68E-02(0.95) 

3.33E-02(0.90) 

c2h3 

1.08E-03(0.80) 

9.05E-04(0.85) 

4.55E-04(0.62) 

3.85E-04(0.70) 

6.74E-0.5(0.67) 

c3h3 

6.77E-04(0.90) 

4.58E-04(0.90) 

2.82E-03(0.90) 

1.01E-03(1.00) 

4.26E-04(1.03) 

c3h5 

6.00E-04(0.70) 

2.88E-04(0.75) 

2.16E-05(0.50) 

6.57E-04(0.65) 

6.80E-0.5(0.78) 

c4h3 

1.23E-04(0.98) 

1.36E-04(1.12) 

6.84E-05(0.90) 

9.85E-05(0.95) 

6.70E-06(0.71) 

c4h5 

2.42E-0.5(0.62) 

8.93E-05(0.70) 

2.79E-05(0.55) 

3.18E-04(0.70) 

2.90E-06(0.58) 

c6h5 

2.08E-06(0.72) 

1.96E-06(0.80) 

3.92E-07(0.06) 

1.62E-06(0.75) 

1.70E-0.5(0.81) 

c6h6 

9.84E-06(0.70) 

1.90E-05(0.45) 

3.54E-06(0.45) 

1.18E-04(0.65) 

3.33E-0.5(0.93) 

(c)  C3H6 

flame — F3 

Species 

Ml 

Mia 

M2 

M3 

Data 

c2h2 

3.29E-02(0.45) 

5.08E-02(0.45) 

5.22E-02(0.56) 

5.79E-02(0.39) 

3.69E-02(0.43) 

c2h3 

6.70E-04(0.35) 

1.10E-03(0.38) 

2.4.5E-04(0.44) 

3.01E-04(0.35) 

2.48E-03(0.34) 

c3h3 

3.95E-03(0.38) 

1.83E-03(0.38) 

7.77E-03(0.47) 

2.88E-03(0.41) 

3.09E-03(0.34) 

c3h5 

1.09E-02(0.31) 

6.07E -03(0.28) 

3.63E-03(0.22) 

1.02E-02(0.32) 

3.50E-03(0.29) 

c4h3 

6.54E-04(0.45) 

7.40E-04(0.40) 

2.78E-04(0.48) 

3.70E-04(0.45) 

1.55E-04(0.43) 

c4h5 

1.01E-04(0.35) 

6. 12E -04(0.29) 

1.04E-06(0.42) 

1.56E-03(0.39) 

1.20E-04(0.34) 

c6h5 

1.06E-04(0.35) 

4.03E-05(0.35) 

1.31E-07(0.41) 

2.86E-05(0.38) 

7.43E-0.5(0.29) 

c6h6 

6.23E-04(0.32) 

1.38E-03(0.25) 

1.08E-05(0.19) 

3.22E-03(0.31) 

1.22E-03(0.39) 

Note:  Numbers  in  parentheses  give  the  peak  location  in  cm  height  above  burner.  Modeling  results  for  many  species 
are  sums  of  predicted  concentrations  of  more  than  one  isomer,  “xxx”  means  that  no  data  were  reported  for  that  species. 
For  FI,  C2H2  is  the  fuel;  the  only  values  reported  are  the  heights  above  burner  at  which  the  mole  fraction  is  half  of  its 
inlet  value. 


observation  of  a  C3H3  +  CH3  product  channel  [30], 
the  branching  ratio  was  set  to  10%  i-C4H5  +  H,  20% 
C3H3  +  CH.3,  70%  C2H2  +  C2H4.  Rates  for  R8 
were  from  Tsang  and  Hampson s  [31]  20  Torr  results. 
C3H  chemistry  was  added,  with  thermochemical 
properties  from  Guadagnini  et  al.  [32],  Allyl/propar- 
gyl  recombination  (R4)  in  M3  formed  much  of  the 
benzene  in  F2  and  F3,  so  this  reaction  was  added, 
although  the  authors  [9]  stated  the  rate  coefficient 
was  an  upper  bound  at  1  atm.  Also  added  were  two 
reactions  for  converting  fulvene  into  benzene.  (See 
Table  1  for  rate  coefficients  of  key  reactions.) 


Radical-Radical  Recombination  Reactions  for 
C6H6  Formation 

Several  chemically  activated  radical-radical  com¬ 
binations  were  proposed  as  possible  benzene-form¬ 
ing  reactions  in  Mia  (R5,  R6,  R9-R12).  For  pro- 
pargyl-propargyl  recombination,  not  only  the  initial 
adduct  but  also  several  of  the  rearrangement  prod¬ 
ucts/intermediates  are  aliphatic  C6H6  species 
[2,5,6];  many  other  radical-radical  combinations  can 
lead  to  the  same  aliphatic  C6H6  species.  Both  n- 
C4H3  +  C2H3  and  n-C4H5  +  C2H  can  combine  to 
form  CH2CHCHCHCCH.  CH3CHCCHCCH, 
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TABLE  3 

Standard  heats  of  formation  (in  kcal/mol)  of  C6H6 
species  and  selected  combinations  of  molecules  with 
combined  empirical  formula  of  C6H6 


Molecule(s) 

AH°f  (298  K) 

n- C4H5  +  C2H 

218.7 

i-C4H5  +  C2H 

206.3 

o-benzyne  +  2H 

205.8 

U-C4  H3  +  C2H3 

201.0 

i-C4H3  +  C2H3 

182.3 

C3H3  +  C3H3 

172.0  [33,34] 

n-C5H3  +  CH3 

169.9 

*-c5h3  +  ch3 

163.1 

c5h2  +  ch4 

147.3 

phenyl  +  H 

131.5 

C4H2  +  c2h4 

123.8 

C4H4  +  C2H2 

123.3 

HCCCH2CHCCH2 

104.3 

CH3CHCCHCCH 

99.7 

CH2CCHCHCCH2 

96.5 

CH3CHCCCCH2 

93.2  [11] 

CH2CCCHCHCH2 

85.6  [11] 

CH2C(CCH)CHCH2 

82.4 

ch2chchchcch 

80.0 

Fulvene 

52.5 

Benzene 

19.8 

Thermodynamic  properties  are  from  Miller  and  Melius 
[5]  and  Kee  et  al.  [20]  unless  otherwise  noted. 


formed  by  H-C5H3  +  CH3,  which  can  undergo  a 
hydrogen  shift  to  form  CH2CHCHCHCCH.  Com¬ 
bination  of  i-C4II  j  +  C2H3  or  i-C4H5  +  C2H  can 
form  CH2C(CCH)CHCH2.  i-C4H5  +  C2H  can  also 
form  CH2CCHCH2CCH,  as  can  C3H3  +  C3H3. 
CH3CHCCCCH2  from  /-C5H3  +  CH3  can  undergo 
a  hydrogen  shift  to  form  CH2CCCHCHCH2. 
CH2CCCHCHCH2  is  also  a  product  of  i-C4H3  + 
C2H3,  converting  via  a  hydrogen  shift  to 
CH2CCHCHCCH2,  another  direct  propargyl-pro- 
pargyl  recombination  product.  Once  any  of  these  al¬ 
iphatic  C6H6  adducts  is  formed,  the  chemistry  lead¬ 
ing  to  either  fulvene,  benzene,  or  phenyl  +  H  is  as 
described  by  Miller  and  coworkers  [2,5,6],  The  over¬ 
all  rate  for  adduct  formation  and  subsequent  reac¬ 
tion  for  each  of  these  new  reactions  was  estimated 
to  be  2  X  1013  cm3/mol  s.  H-abstraction  reactions 
from  these  pairs  of  radicals  are  also  included  in  the 
new  mechanism,  as  are  reactions  interrelating  CsH2, 
i-C5H3,  and  n-C5H3. 

All  these  new  reactions  were  written  in  the  exo¬ 
thermic  direction.  Standard  heats  of  formation  for 
the  above  C6H6  species  and  combinations  of  mole¬ 
cules  with  combined  empirical  formula  of  C6H6  are 
given  in  Table  3.  The  heat  of  formation  used  for 
propargyl  radical  is  86.0  kcal/mol  [33,34],  and  the 


standard  entropy  is  59.9  cal/mol  K  [34];  the  earlier 
values  used  were  83.0  kcal/mol  and  61.4  cal/mol  K. 
The  less  stable  propargyl  radical  results  in  slightly 
lower  predicted  C3H3  mole  fractions,  and  therefore 
a  slightly  diminished  role  of  C3H3  +  C3H3  in  ben¬ 
zene  formation.  Predicted  peak  mole  fractions  of 
C4H3  and  C4H5  also  increased  by  at  least  a  factor  of 
3  in  all  three  flames.  The  total  rate  used  for  R3  was 
1  X  1013  cm3/mol  s;  measured  values  were  2.4-7.2 
X  1013  cm3/mol  s  [35-38].  Preliminary  master 
equation  calculations  [39]  showed  that  at  conditions 
of  the  flames  considered  here,  veiy  little  benzene 
was  formed  directly.  Accordingly,  the  branching  ratio 
used  had  the  products  being  split  nearly  evenly  be¬ 
tween  fulvene  and  phenyl  +  II,  with  i-C5H3  +  CH3 
being  a  minor  pathway  since  it  is  slightly  exothermic. 
(Not  enough  confidence  was  placed  in  the  2.1  kcal/ 
mol  exothermicity  of  the  n-C5H3  +  CPI3  product 
channel  to  believe  that  it  would  be  significant.) 

Discussion  of  Results  from  Improved 
Mechanism 

Predictions  using  Ml  and  Mia  are  compared 
against  data  in  Fig.  1  for  CgHe.  Benzene  is  the  only 
CgH6  species  in  Ml;  Mia  also  includes  fulvene.  For 
FI  (Fig.  la),  Mia  greatly  improves  C6H6  predic¬ 
tions.  Little  fulvene  is  predicted.  An  improvement  is 
also  seen  for  F2  (Fig.  lb),  with  an  appreciable 
amount  of  fulvene  predicted  early  in  the  flame.  Mia 
slightly  overpredicts  C6H6  concentrations  in  F3  (Fig. 
lc),  with  fulvene  being  important  early  in  the  flame. 

The  major  reaction  pathways  leading  to  phenyl, 
fulvene,  and  benzene  are  shown  in  Fig.  2.  For  all 
three  flames,  R3  plays  a  major  role.  R4  contributes 
a  moderate  amount  in  F2  and  F.3;  in  FI,  it  is  a  major 
fulvene  decomposition  pathway.  Use  of  the  1  atm 
rate  coefficient  for  this  reaction  is  most  likely  re¬ 
sponsible  for  the  overprediction  of  C6H6  in  F3.  R1 
is  almost  entirely  a  phenyl  decomposition  pathway 
in  F2  and  F.3;  in  FI  the  reaction  is  of  negligible 
importance.  R2  has  no  discernible  effect  for  any  of 
the  flames  studied. 

Of  the  new  radical-radical  reactions  added,  only 
R5a  appears  as  a  substantial  ring-forming  reaction, 
and  then  only  in  FI.  The  relation  of  the  C5FL  and 
C5H3  species  with  the  polyacetylene  pool  leads  to 
more  of  these  species  being  formed  in  FI.  Unfor¬ 
tunately,  there  are  no  data  for  testing  predicted  con¬ 
centrations  of  the  CsH2  and  C5H3  species.  For  F2 
and  F3,  the  major  role  of  R5  is  as  a  pathway  for 
converting  fulvene  to  benzene.  Whether  this  actually 
happens  is  questionable;  however,  further  investi¬ 
gation  of  the  possible  role  of  these  CsH3  species  is 
indicated.  The  only  other  new  reactions  appearing 
in  Fig.  2  involve  j-C4H3  +  C2H3  (Rll)  as  decom¬ 
position  products  of  either  fulvene  or  phenyl  +  U. 

No  C2  +  C4  reaction  in  the  mechanism  was  pre¬ 
dicted  to  lead  to  benzene  formation  in  any  of  these 
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Fig.  1.  Predictions  of  Cfl 1 1 g  species  from  mechanisms 
Ml  and  Mia  compared  with  data:  (a)  C2H2  flame  (FI),  (b) 
C2H4 flame  (F2),  (c)  C3H6  flame  (F3).  Circles,  data;  dashed 
line,  benzene  from  Ml;  open  squares,  benzene  from  Mia; 
open  triangles,  fulvene  from  Mia;  solid  line,  total  C6H6 
from  Mia. 


flames — only  reactions  involving  C3  radicals  led  to 
ring  formation.  Fig.  3  shows  the  major  predicted 
pathways  for  C3  radical  formation  from  the  fuel  for 
the  three  flames.  For  F3,  the  main  source  of  C3H3 
and  C3Fl5  radicals  is  by  successive  H  abstractions 
from  the  fuel.  While  R7  is  the  dominant  source  of 
C3  radicals  in  FI  and  F2,  it  is  a  secondary  pathway 
in  F3,  with  the  fuel  decomposing  to  ethylene  or  vinyl 
radical,  forming  acetylene  by  H  abstraction.  In  F2, 
C2H4  and  C2H3  combine  with  C4  radicals  to  form 
allyl.  In  Fl,  the  fuel  combines  with  methyl  radical 
to  give  propyne  +  H,  with  the  propyne  then  forming 
propargyl. 

Therefore,  for  the  low-pressure  flames  studied 
here,  only  combinations  of  C3  radicals  led  to  pro¬ 
duction  of  benzene  or  phenyl.  However,  it  is  con¬ 
ceivable  that  there  would  be  fuels  or  flame  condi¬ 
tions  that  would  favor  the  C2  +  C4  pathways.  For 
example,  R2  is  the  only  important  molecular- weight 
growth  reaction  involving  phenyl  having  a  net  mole 
change.  The  rates  of  bimolecular  reactions,  ex¬ 
pressed  as  mol  C6H5/g  gas/s  (the  molar  rate  divided 
by  the  gas  density),  would  increase  linearly  with 
pressure,  while  those  of  unimolecular  rates  would 
remain  constant,  all  else  being  equal.  The  relative 
importance  of  R2  to  the  other  molecular-weight 
growth  reactions  at  higher  pressures  can  then  be  ap¬ 
proximated  by  looking  only  at  its  forward  rate,  which 
is  2.4%  of  the  peak  of  the  sum  of  the  net  rates  of  all 
the  other  molecular-weight  growth  reactions  for  Fl, 
and  less  than  1%  for  F2  and  F3.  Therefore,  in¬ 
creased  pressure  alone  would  decrease  the  role  of 
R2  as  a  phenyl  decomposition  channel  but  would  not 
in  itself  make  R2  a  major  phenyl-producing  reaction. 
This  is  more  noteworthy  since  C4H3  was  overpre¬ 
dicted  in  all  three  flames,  and  since  the  n  —  C4H3/ 
i  =  C4H3  and  n-C4H5/i-C4Hs  ratios  are  at  least  an 
order  of  magnitude  above  partial  equilibrium  for  all 
positions  in  all  three  flames.  Secondary  effects  of 
pressure,  such  as  increasing  the  rates  of  C4-radical 
forming  reactions,  cannot  be  ruled  out. 

In  summary,  for  the  three  flames  considered,  the 
only  major  formation  pathways  for  benzene  or 
phenyl  were  the  result  of  combination  of  C3  radicals. 
The  C2H2  +  C4  radical  pathways  (Rl,  R2)  did  not 
appear  to  contribute  significantly  to  single-ring  aro¬ 
matic  formation.  The  relative  importance  of  differ¬ 
ent  benzene/phenyl/fulvene  reaction  pathways  was 
fuel  dependent,  as  were  the  major  pathways  for  for¬ 
mation  of  the  C3  radicals.  No  claims  to  generality 
can  be  made  for  other  fuels,  flow  systems,  or  pres¬ 
sures.  However,  in  modeling  of  low-pressure  C2H2 
flames  [40]  and  atmospheric-pressure  C2H4  flames 
[40-42],  recombination  of  propargyl  radical  was 
found  to  be  an  important  contributor  to  benzene  for¬ 
mation. 

The  new  pathways  proposed  to  C6H6  are  pertinent 
to  the  overall  role  of  RSFRs  in  benzene  formation, 
since  the  aliphatic  C6H6  adducts  formed  are  the 
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H2CCCH 


Fig.  2.  Predicted  major  reaction  pathways  involving  phenyl,  benzene,  and  fulvene  for  the  three  flames  from  mechanism 
Mia.  See  text  for  discussion. 


Fig.  3.  Predicted  primary  and  sec¬ 
ondary  formation  pathways  of  C3 
radicals  for  the  three  flames  from 
mechanism  Mia.  See  text  for  discus¬ 
sion. 
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same  as  for  the  C3H3  +  C3H3  reaction.  The  C5H3 
+  CH3  reactions  played  a  role  in  the  predicted  ben¬ 
zene/phenyl  formation  chemistry,  perhaps  due  to 
their  being  only  ~32-38  kcal/mol  higher  in  energy 
than  phenyl  +  H.  However,  due  to  concentration 
effects  and  a  net  increase  in  number  of  moles,  even 
the  C2  1 1 2  +  /-C4H3  channel,  being  50.8  kcal/mol 
higher  in  energy  than  phenyl  +  H,  was  a  phenyl 
decomposition  channel  in  FI.  There  is  a  need  for 
clarification  of  the  actual  rates  and  branching  ratios 
of  these  unsaturated  radical-radical  combinations 
that  can  undergo  multiple  rearrangements,  espe¬ 
cially  considering  some  of  the  possible  new  pathways 
presented  here. 

Consideration  of  formation  pathways  of  larger  ar¬ 
omatic  species  can  be  expected  to  share  the  same 
type  of  effects  of  the  fuel  and  flame  conditions.  For 
example,  in  the  present  work,  essentially  all  of  the 
naphthalene  formed  in  F2  and  F3  was  predicted 
(from  Mia)  to  be  formed  by  C5H5  +  C5H5.  In  FI, 
this  reaction  was  second  in  importance  to  C8H5(o- 
ethynylphenyl)  +  C2H2  <=>  C10H7  (1-naphthyl).  For¬ 
mation  of  larger  PAHs  can  reasonably  be  expected 
to  be  less  strongly  affected  by  falloff  behavior,  due 
to  the  larger  molecules’  ability  to  absorb  the  energy 
released  by  chemical  activation.  It  looks  as  if  the 
most  complete  answers  to  questions  of  the  formation 
of  aromatic  compounds  might  never  be  simple  ones. 
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COMMENTS 


K  -ll  Homann,  TU  Darmstadt,  Germany.  The  very  early 
maximum  concentration  of  C2H,  both  in  acetylene  and  eth¬ 
ylene  flames,  could  be  somehow  connected  to  the  forma¬ 
tion  of  benzene  that  is  also  a  rather  early  product.  Have 
you  any  new  ideas  how  C2H  could  be  formed  so  early,  that 
is,  in  the  lower-temperature  part  of  the  oxidation  zone?  The 
traditional  path  H  +  C2H2  — >  C2H  +  H2  is  inadequate  to 
explain  the  early  maximum  of  C2H. 

Author’s  Reply.  In  the  acetylene  flame  considered  here 
(Ref.  [16]  in  paper);  neither  benzene  nor  C2H  can  really 
be  considered  an  early  product,  with  the  peak  measured 
mole  fractions  occurring  well  into  the  reaction  zone  (at  0.6 
and  0.77  cm,  respectively),  and  no  significant  amount  of 
either  before  roughly  0.3  and  0.4  cm,  respectively.  Mia 
predicts  the  peak  locations  well  (0.625  and  0.9  cm,  respec¬ 
tively).  Benzene  plays  no  apparent  role  (direct  or  indirect) 
in  C2H  formation.  The  main  formation  reaction  is  C2H2  + 
OH  C2H  +  H20,  which  is  roughly  5  times  more  im¬ 
portant  than  C2H2  +  II  < — >  C2H  +  H2.  Alternately,  in  the 
ethylene  flame  (Ref.  [17]  in  paper);  the  measured  C2H 
peak  occurs  at  0.83  cm,  close  to  the  benzene  peak  (0.90 
cm).  While  mechanism  Mia  predicts  well  the  measured 
C2H4,  02,  and  C2H2  profiles,  not  only  does  it  predict  too 
early  a  peak  benzene  mole  fraction  (0.45  cm),  the  peak 
predicted  C2H  mole  fraction  is  too  late  (1.5  cm).  Again, 


CoH2  +  OH  <-»  C2H  +  H20  is  the  major  predicted  C2H 
formation  reaction,  but  C2H2  +  H  C2H  +  H2  becomes 
a  net  C2H  destruction  reaction.  The  C4  reactions  described 
above  lead  to  net  production  of  C2H,  becoming  more  im¬ 
portant  later  in  the  flame.  The  possible  discrepancies  be¬ 
tween  the  data  and  model  for  C2H  might  therefore  reflect 
limitations  in  the  C4H3  chemistry.  There  is  also  a  predicted 
contribution  to  C2H  formation  from  C3H3  t  ()o  CH20 
+  C2H.  Since  both  C3H3  and  G4H3  are  relevant  to  ben¬ 
zene  formation,  possible  significant  coupling  of  the  ben¬ 
zene  and  C2H  chemistry  cannot  be  ruled  out. 


Katharina  Kohse-Hoinghaus,  University  of  Bielefeld, 
Germany.  Does  your  updated  model  include  NO,-related 
chemistry  under  fuel-rich  conditions,  or  are  you  consider¬ 
ing  to  include  this?  It  would  be  very  helpful  if  NO,  re¬ 
burning  chemistry  for  other  fuels  than  methane  could  be 
reliably  modeled. 

Author’s  Reply.  J.  A.  Miller  has  developed  a  reaction  set 
for  describing  the  C/H/O/N  chemistry,  which  is  available 
from  the  authors  upon  request.  No  changes  in  the  NO, 
chemistry  resulted  from  the  present  work,  since  this  was 
not  a  focus  of  the  work,  and  none  of  the  flames  studied  in 
the  present  work  even  contain  nitrogen. 


